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I. INTRODUCTION 

Precise measurement of CP violating parameters in weak decays is an active topic. On the other hand, testing 
CPT symmetry is also of great importance, as some physics beyond the standard model may lead to its breaking. 
Nowadays, neutral K and B mesons are routinely produced as quantum-entangled or EPR-correlated pairs in 4> and 
B factories (H-H)]. Hence it is highly interesting to explore their use in examining CP and CPT symmetries [JH, [H~ 
[To| . and in measuring relevant parameters (ill . H2| . Especially, systematic proposals were made on exploring CP 
and CPT violating parameters in correlated decays of meson pairs in C — — 1 state, by using decay rates and their 
asymmetries @, 0|- In this paper, we present a general and rigorous calculation of time-ordered integrated rates of 
the correlated decays of C = ±1 entangled meson pairs, taking into account both CP violation and possible CPT 
violation. Subsequently CP and CPT violating parameters are expressed as functions of the four asymmetries defined 
for the time-ordered integrated rates of the correlated decays of the two entangled states. Our result provides a general 
method of measuring rephase- invariant CP and CPT violating parameters in correlated decays of entangled mesons 



into any kind of product. This extends a rephase-invariant formalism of CP and CPT violating observables 13[ to 
C = ±1 entangled mesons. 

In Sec. |nl we review the single-meson decays using two complex numbers parameterizing indirect CP and CPT 
violations. In Sec. lIIIl we consider the entangled states \i/j±} and calculate the general expressions of the decay rates as 
functions of four single-meson decay amplitudes. In Sec. IIV1 we integrate the decay rates over all times with two time 
orders, thereby obtain four time-ordered integrated rates. We give both exact results and approximate expressions 
up to the first order of the CP and CPT violating parameters. Subsequently, we obtaining the four real parameters 
as functions of the four asymmetries defined for the four time-ordered integrated rates. In Sec. [Vj two special cases 
are treated for which the expressions are much simpler than the general expressions in Sec. IIVI Discussions and a 
summary are made in Sec. I VI I 

II. REVIEW OF CP AND CPT VIOLATING PARAMETERS IN SINGLE-MESON DECAYS 

An arbitrary state of a pseudoscalar meson can be written as 

\M{t))=m\M°)+m\M°), a) 

where \M°) is the pseudoscalar meson state and |M°) is its antiparticle state, ip(t) and -ip(t) are superposition 
coefficients. In Wigner-Weisskopf approximation, the time evolution of \M(t)) is described as 



.±fm\ = f h u h 12 \ fm\ 

% dt\${t)) \H 21 H 22 J \i>(t)J> 
where Hij — Mij — iTij/2, with the dispersive part and the absorptive part are hermitian. Define [lJ 



(2) 



H 22 - Hu 

6 M = j= TT , (3) 



\JH\ 2 H 2 



q _ H 21 _ l - e M 



The eigenstates of the Hamiltonian are 



p V H 12 1 + e M 

\M 1 )=p 1 \M°) + q 1 \M ), 
\M 2 ) =p 2 \M°) -q 2 \M°), 



(4) 



(5) 



with the respective eigenvalues 



Ai — Hu + y H\ 2 H 2 i{ 







S M 




4 


Sm , 




SI, 

4 



(6) 



Therefore we obtain 



| M o } = Q2\Mi) + qi\M 2 ) 



\M°) 



viqi +P2qi 

p 2 \Nh) -pi\M 2 ) 

pi 92 +£>2<2 , l 



(7) 
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For a single meson whose initial state is |M°), we have 

mt)) = Plg2e " Alt+P2gie " a2 Vo) + ^ e ' iXlt - e ~ iX2t) \M ). (8) 
Piq2+P2qi piq2+P2qi 

Similarly, for a single meson whose initial state is |M°), 

mt)) - piP2{e ~ iXlt - e ~ iX2t) \M 0) + p^- iXit +p^- iX2t m . (9 ) 

Piq2+P2qi piq2+P2qi 

Define 



1 + A M _ S M . L , S 2 M 



1 - A M 2 



then hi 



Pi pi- A 



M 



qi ql + A M ' 
P2 _ P 1 + A M 
g 2 5 1 - Am ' 



Therefore 



V2 
1 

i 

Therefore, up to the order of O(Am) and 0{em), we can simplify Equations (JSJ and (|9|) as 



where 



(10) 

(11) 
(12) 



pi oc -L(i+ £M )(i- A M ), (13) 
<ji oc — (l-e M )(l + A M ), (14) 



;m -v — =(1 - - A.w). (!•')) 

92 cx _(i- eM )(i- A M ). (16) 



\M°(t)) = [g + (t) - 2A M g_(t)]\M ) + (1 + 2e M )g_(t)\M ), (17) 
\M°(t)) = (1 - 2e M ).9- (t)|M ) + [<?+(*) + 2A M ff-(*)]|Mo}, (18) 



.g + ^i( e - lAlt + e -^ 4 ), (19) 
ff _ = i( e -^*- e -^*). (20) 

If CPT is conserved, then Hu — H22, therefore Am = 0. Independent of the situation of CPT, if CP is conserved, 
then £m =0. If T is conserved, then M^ 2 /M 12 — r* 2 /r 12 , therefore 5?£m = 0. If CPT and CP are both conserved, 
we have T conservation. But T conservation can also be valid without requiring the conservation of CPT and CP. 
T conservation together with CP violation means 3sm 7^ 0. 

III. CORRELATED DECAYS OF ENTANGLED STATES 

Suppose the initial state of two mesons a and b is the entangled state of C = ±1, 

|V±> = -j=[\Mo) a \M ) b ± \M ) a \Mo) b }. (21) 



After time t a , the final state of a is |/ a ), while after time %, the final state of b is \fi 



b ■ 
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Define the decay amplitudes of Mq and Mq 



91 = (f a \M ), (22) 

92 = (fb\M ), (23) 
hi = (f a \M ), (24) 
h 2 = (f b \M ). (25) 



The decay amplitude of the entangled meson pair can be obtained from ([T71) and ([T8|) as 



A lh (*«'**) = ^[(/-l M o(ia))(/ 6 |Mo(4)) ± (/«|Mo(ta))</6|Af (t6))]. (26) 



The decay rates for the entangled states are thus 
Tl h (t a ,t b ) = \A+ afb (t a ,t b )\ 2 



where 



i e -rx(*.+t fc )| 0l |2 + I e -r a (t.+t fc )| Q2 |2 + l e -r lto -r 2tb|a3|2 + l e -r,t a -r lt > ^2 

+e- r( '' +il »SKa 2 e- ,Am(t « +t,) ] + e- r ^- rtb ^[ala 3 e- lAmtb ] 
+e- rt «- Tlt '>%i[a* 1 a 4 e- lAmta } + e- rt *- V2tb $l[a* 2 a 3 e iAmt «} 

+e- r2ta - rtb R[a* 2 a A e iAmtb ] + e^^+^^a^e^" 1 ^^], (27) 



r 



+e- r ^ +tb ^lp2e- iAm ^ a+tb) } + e- T ^- Ttb ^.[l3l^e- lAmtb ] 
+e- rt "- rit ^[/3*/3 4 e- lAmt »] + e- Tt °- T ^R[(5* 2 foe tAmt "] 

+ e - r2t »- rt ^[/3 2 */3 4 e lAmtb ] + e- v{ - ta+tb ^[^ A e iAmi - tb - ta \ (28) 



Am = m 2 — mi, 

r = Ei±Ii 

2 ' 

Ar = r 2 -ri, 

on = 2(51 + h i)(92 + h 2 ) + (A M + e M ){hih 2 - gig 2 ) - (A M + e 2 M )(gih 2 + g 2 hi) 

+2A M e M (gig2 + hih 2 ), 
a 2 = -\(gi - h\){g 2 - h 2 ) + (A M ~ e M )[hih 2 - gig 2 ) - (A M + £ 2 M )(gih 2 + g 2 hi) 

+2A M e M (g 1 g 2 +h 1 h 2 ), 
a 3 = A M (.91 + hi)(g 2 - h 2 ) + (A 2 M + e 2 M )(gih 2 + g 2 hi) - 2A M e M (gig2 + hih 2 ), 
a 4 = A M (gi - hi)(g 2 + h 2 ) + (A 2 M + e 2 M ){gih 2 + g 2 hi) - 2A M e M (gig2 + hih 2 ), 
p 1 =fo = (A M - s M )(g 2 h! -gih 2 ), 

fa = 2(51 + hi)(h 2 - g 2 ) - A M {g 2 h\ + g\h 2 ) + e M {g\g2 + h\h 2 ) + (A 2 M - e 2 M )(gih 2 - g 2 h x ), 
At = 7}(gi - hi)(g 2 + h 2 ) + A M (g 2 hi + g\h 2 ) - s M (gig2 + hih 2 ) + (A 2 M - e 2 M ){gih 2 - g 2 hi). 

IV. TIME-ORDERED INTEGRATED DECAY RATES 

Now consider two kinds of decay of an entangled state \tp±). One, whose integrated rate is denoted as Rf a , is that 
the decay of b into fb precedes that of a into f a , i.e. t b <t a . The other, whose integrated rate is denoted as R^ b , is 



of the inverse time ordering, i.e. tb> t a . That is, 



poo pt a 

R ba = dt a dt b T± (t a ,t b ), 

Jo Jo 

/>oo poo 

R ab = dt a dt b Tf f (t a ,t b ). 

JO Jt„ 



Making use of Equations ([27)) and (|28|) . we obtain 

jy+ l«i| 2 , M 2 , |«3| 2 , M 2 . aja^ 1 , ^ a\a 3 

R ba - TFT + TFT + TFT + TFT + ^Ufr , „-a^ 2 J + K Li 



4rf 4r| 4rir 4r 2 r L 2(r + iAm) 2J L r x (ri + r + iAm) 



(r + iAm) (Ti + T + iAm) J L (r - iAm) (T 2 + T - iAm) J 

^4 n R r 0|04 n 

L r 2 (r 2 + r-zAm) J L 2r(r + iAm) J 

~ A + Ai$RA A/ + A 2 SA M - A^te M - A£se M , 

n . |/?l| 2 . |fe| 2 . |fe| 2 . |/?4| 2 1 r ftfh 

n ba Ar 2 Ar 2 /1F.F /IF„F Lo/T' _L h'AotoAS J Ll 



4r 2 4r| 4rir 4r 2 r L 2(r + iAm) 2J L ri(ri + r + iAm) J 



(r + iAm)(ri + r + iAm) J L (r - iAm)(r 2 + T - iAm)- 



-T 2 (r 2 +T -iAm) J L 2r(r + iAm) J 
B + BxSRAm + B 2 3A M - B 3 $te M - B 4 $e M , 

l^ll 2 , K| 2 |Q3| 2 |Q4| 2 Qj«2 i s^r Olf_3 i 

4r 2 4r 2 + 4r 2 r 4rir L 2(r + iAm) 2j L (r + iAm)(ri + r + iAm) J 



a 1 a 4 -1 r a 2 o;3 



r x (ri + r + iAm) J L r 2 (r 2 + r-iAm) J 

a 2 a 4 1 crvf 0^30:4 



+SR[ _^4 1 

L ( r -* Am )( r 2 + r-iAm) J L 2r(r-zAm) J 
~ Co + CxKAm + C 2 9A M - C 3 Ke M - C 4 ^ M , 

d- _ 1 IfeP , Ifel 2 . I&l 2 . ^ flft 1 

■"-ab /1F2 ^r2 /IF„F /IF.F L Ol'F x_ „' A ™ "|2 J L / 



4r 2 4r 2 4r 2 r 4rir ' L 2(r + iAm) 2J L (r + iAm)(ri + r + iAm)- 



r x (Ti + r + iAm) J L r 2 (r 2 + r-iAm) J 



■(r-iAm)(r 2 + r-iAm) J L 2r(r~iAm) J 
Co + £>iKA M + ^2^A M - £> 3 3? £M - D^sm, 



where {Ai}, {Bi}, {d} and {Di}(i — 0, 4) are given in the appendix. 

Now we consider the asymmetries defined for these four time-ordered integrated decay rates, 



R ba 


- R lb 


R ba 


+ R tb 


R ba 


- R ab 


R ba 


+ K b 


R ba 


- R ba 


R ba 


+ R ba 


R th 


— R nh 
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which are calculated to be 

Ri 
R2 



M - %A M + M 



C 2 



2A a 2A 
B -D Q 2(B 1 D Q -B Q D 1 
Bp + D (Bo + D ) 2 



3 Am, 



2(B 2 Dq Bp D 2 ) 

■^ Am+ {b + DoY " Am 



H — — ^ ^ »i£m H 7tz — , ,-. s£m 1 



(Bo + Dof 



(B + Do) 2 



Aq-Bq 2(A 1 Bq-A B 1 ) ^ a 2(A 2 Bo-AqB 2 ) ^^ 
R3 — — — — — H — — „ „ 9iAm H — — „ xo <sAm 



+ Bo (Ao + Bo) 2 



(A Q + Bo) 2 



2(A 3 B -ApB 3 ) 2(A 4 B - A B 4 ) ^ 
H — — „ NO Uem H — — , ^9 — s£m, 



(Ao 



(Ao + Bo) 2 



Ha 



which can be rewritten as 



Co 



Do 2(dA) -CqD 1 ) soa , 2(C 2 Dq-CqD 2 ) ^ a 
H — — „ ^ St^Af H 77=; — , „ N9 — 



Co + £ (Co + A)) 2 



(Co + Do) 2 



2(C 3 D - CqD 3 ) ^ 2(C 4 D - C D 4 ) ^ 



(Co + Do) 2 



(Co + D ) 2 



( 



Ri 

_ Bo --Do 
B0+-D0 



i? 



3 A +B 



K 



with 



/ 



K = 



2An 



-Ci 



-c 2 



Z (A +B ) 2 



2A 

q -BiJo-BqDi r, B 2 D -B D 2 o B 3 D -B D 3 r, B 4 D„-B„D 4 
Z (B +B ) 2 Z (B + B ) 2 - ^ -.,>,,,,■> 







o A 2 -B -A B 2 
Z (A +B ) 2 





- »o- .. 
(Bo+Bo) 2 



O A 4 j?Q— A j? 4 
Z (Ao+R^ 2 



o Ci-Dp — CqD\ q C2-D0 — CqD2 q C3D0— C0D3 o C4D0— C0-D4 
\ Z (C0+B0) 2 Z (C0+B0) 2 Z (C0+B0) 2 Z (C +B ) 2 / 

Therefore, CP and CPT violating parameters are obtained as 

\ 



/ M M \ 








V 3e M / 





/ 

R 2 
R 3 



Ri 

_ Bp- D 
Ba+Da 

_ Aq-Bq 

Co+Bo 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



where, in terms of the Levi-Civita symbol tijku ^-ijk = (AiB — A Bi)(BjDo — B Dj)(CkD — CoDk) and Ayk 
(A 2 — C 2 )e i2 jk — (Ax — Ci)eujk(i,j, k = 1, 2, 3, 4), the matrix elements of K~ x are 



{K-% 2 
(K-% 2 



2Ao€iijk^ijk 

^ijk^-ijk 

(A 2 - C 2 )(A 324 - A 423 )(B + Dp) 2 

2(B 2 Do — BoD 2 )\ijtAijk 
(A 2 - C 2 )(A 234 - A 243 )(A + Bp) 2 
2(A 2 B — A B 2 )\ijkAijk 
(A 2 - C 2 )(A 432 - A 342 )(Cp + Dp) 2 
2(C 2 D — CpD 2 )XijkAijk 
2ApEi 2 jkAijk 
X A ' 

(Ax - Ci)(A 3 i4 - A 413 )(Bp + Dp) 2 
2(B\Do — BpD\)\ijkAijk 
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(K- 1 
(K- 1 



23 — 



24 



31 



32 — 



43 



44 



(A 1 -C 1 )(A 134 -A 143 )(A + B o y 
2(^4i£?o — AoBi)XijkAijk 

(A 1 -C 1 )(A 4 3l-A34l)(Co + ^o) 2 

2(CiD — Co-DijAyfeAyfc 



AjjfcA 



ijk 



33 



34 — 



11 



42 — 



17 A 

[(A432 


A \( A 
— A 2 34)(Ai 


- Oi) - 


/A A \ / A 




n \2 
















A \! A 


- Ci) - 


1 A A \ / /l 
(A 3 i 4 - A 3 4iJ(A 2 


- C 2 )\\Aq - 


1- -DO J 






^(,^3 -DO 


4 R„ s l\.., A .., 






[7 A 

[(A 2 43 


- A 4 2 3 )(Ai 


- Ci) - 


/A A \ / /I 

(A143 - A 4 i 3 J(A 2 


- <^2)J(C-0 ^ 


n \2 






2(C 3 A> 


— C'o£' 3 )AijfcAjjfe 






2Aoeij 


kiAijk 










Xijk 


Aijk 










[(A 243 


- A 342 )(Ai 


-Ci)- 


(A143 - A 341 )(A 2 


-Cb)](Bo- 








2(fl 4 A) 


— B D4)XijkAijk 






[(A432 


- A 4 2 3 )(A! 


-co- 


(A431 - A413XA2 


-C 2 )](Aq- 


1- -Bo) 2 






2(^4^0 


— Ao-B 4 )AyfcAij7 S 






[(A 3 24 


- A 234 )(Ai 


-Cx)- 


(A314- Ai 34 )(A 2 


-C 2 )](C - 


FA,) 2 






2(C 3 Aj 


— CoD 3 )XijkAijk 







In general, CP and CPT violating parameters can be obtained from the four asymmetries of the time-ordered 
integrated rates of the correlated decays of CP = ±1 entangled states, using the above formulae. But in some 
cases, the four asymmetries may be independent oi CP and CPT violating parameters, as in one of the special cases 
discussed below. 



SPECIAL CASES 



We now consider the following two special cases [13|], for which the calculations are simplified a lot. 

First, we consider the situation that f a = f while /& = /, where / and / are mutual CP conjugates and are 
decay products of M° and M° respectively. For example, for M° = B°, f = D~ Dg , D~ K + ,n~ Dg ,ir~ K + while 
/ = D+D S ,D+K-, n+D s , w+K ~, respectively. 

In this case, the amplitudes satisfy 32 = hi = 0. Up to the order O(Am) and O(em), {oi}, are given as 



1 



ai = a 2 



z9ih 2 , 



a 3 = -on = -A M gih 2 , 
Pi =/3 2 = 0, 

(\ 

{7: + Am) 51/12 



Am) 9ih% 



(46) 

(47) 
(48) 

(49) 
(50) 



Then we obtain A 3:4 = £?i )2 , 3) 4 = C 3 , 4 = A, 2,3,4 = 0. Using the expressions given in the appendix, Eg. (1431) is 
simplified as 



/ 

R2- 
R 3 - 



Ri 



Bo 


-D 






t 


-DO 








--Do 



V 



A1-C1 


A2-C2 


2A 


2A 










A 2 B„ 


(Ao+Ba)'^ 


Z (A +So 











i* 



0/ 





/ SJAm \ 








5Rem 1 


) 


V / 



(51) 



Co+Do / 

In this case, we cannot obtain the CP violating parameter from the four asymmetries, which are independent of 
the former. Moreover, R 2 = Bo+d", an< ^ ^ 4 = Co+go are m d e P en dent of both CP and CPT violating parameters, 
while i?i and i? 3 depend only on the CPT violating parameter, but not on the CP violating parameter. 
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Now we consider the situation that f a = h = f, with / = / being a CP eigenstate. For example, f a = fb = 

7r + 7r~ ,7r°7r°. In this case, the decay amplitude satisfies g\ — g 2 and h\ = h 2 . Hence the {«;}, and {ft} can be 
simplified as 

ai = i(.9i + ft) 2 + (A M + e M ){h\ - g\) - 2(A M + e 2 M ) gi hi + 2A M s M (g 2 + /if), (52) 

"2 = — - h x ) 2 + (A M - e M ){h\ - g\) - 2(A M + e M ). 9 ift + 2A M e M (gf + hj), (53) 

a 3 = a 4 = A M ( ff ? - ft?) + 2(A M + £ M ) 5 ifti - 2A M e M (gl + ft?), (54) 

ft = ft = 0, (55) 

ft = -ft = ~(s 2 - hi) - 2A M <?ift + e M {gl + h 2 ) - 2A M g 1 h l + s M (gf + h\). (56) 



Consequently R ba — R ab = 2\em — ^A/| 2 |ffi| 2 (p^ r 2 +Am 2 ) > wn i cn leads to R2 — exactly, no matter whether 

CP or CPT is violated. 

If CP is conserved, then gi — g 2 — h\ — h 2 , hence a\ = 2g 2 — 2g 2 (Aj^i—eM) 2 , «2 = — "3 = — "4 = — 2# 2 (Am— £m) 2 , 
ft = -ft = 2sr?(e M - Am). Consequently, R+ a = R+ b = + 0((e M - Am) 2 ) • Then up to the order 0{e M ) and 
0(A M ), we have R\ = and R 3 — R^m 1. 



VI. DISCUSSIONS AND SUMMARY 



Earlier, extensive discussions were made on [HlZI- However, most of these discussions concern the time-order 
asymmetry for integrated decay rates with a given time difference, rather than R 2 , which is the asymmetry for the 
integrations over all time differences. Only for the case of_/ a = 7r + 7r~ while fb = tt°tt°, was R 2 calculated, with a 
nonzero value dependent on the CP violating parameter [y, 0|- This is consistent with our results in Sec. IIVI Our 
results are in given in terms of the decay amplitudes defined for \Mq) and \Mo), which are related to the indirect CP 
violating parameter e'. 

In (j) factory and in T(4s) resonance in B factories, the branch ratio of C-even meson pairs is negligibly small 0], 
thus R~^ b and R ba cannot be measured there. However, in T(5s) resonance operated in B factories such as CLEO 
and BELLE, the strong decay products include B S B S , B*B S and B S B*. As B* — > B s j, the final states of B*B S and 
B S B* are B S B S pair in It turns out that the branch fraction of |^_) and are 90% and 10%, respectively 

[HQ]. Therefore, Tf . can both be measured once the correlated pairs are identified. In CPLEAR experiment on 
entangled kaons, the branch ratio between and \ip-) was as considerable as 0.037 [l5j]. Moreover, we hope that 
in future, neutral meson production can start with an initial state of even orbital angular momentum so that pairs in 
1^+) dominate. 

To summarize, we have presented a mo del- independent formalism of extracting rephase-invariant CP and CPT 
violating parameters from four asymmetries of time-ordered integrated rates of correlated decays of C — ±1 entangled 
states of neutral pseudoscalar mesons. Asymmetries of time-ordered integrated rates are used because they are easy 
to be obtained from experimental data, without the necessity of measuring the exact times of decays. We give the 
general formulae applicable to all kinds of decay product. In the special case that M° and M° decay to CP conjugates, 
the asymmetries so defined cannot give information on the CP violating parameter. Then one needs to use other 
observables, such as the time-ordered integrated rates themselves, or the asymmetries defined for specified differences 
of decay times. In future work, we will make more phenomenological studies and a more detailed comparison with 
the previous results [f| 0] ■ 
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Appendix A: Detailed expressions of Ai, Bi, d and Di 



_. \9i + hi\ 2 \g 2 + h 2 \ 2 \gi - h^ 2 ^ - h 2 \ 2 

sr 2 + sr 2 

r 2 - Am 2 

4(r 2 + Am 2 ) 2 [ (k|2 - l^l 2 )^ 2 ' 2 ~ ^ - ^(9iKMg 2 h* 2 )] , 

1 = 2f5 [(M^l 2 - I51.92I 2 ) - (|.gi| 2 - IhtfWfahZ) - (\g 2 \ 2 \h 2 \ 2 W( 9l hl)] 
- ^ [l^i^l 2 - \gig 2 \ 2 + (| 5 i| 2 - \hi\ 2 W(g 2 h* 2 ) + (\g 2 \ 2 - \h 2 \ 2 )M{ gi hlj\ 



2 

T 2 - Am 2 



(r 2 + Am 2 ) 



[(M 2 - |/n| 2 M 52 ^) + (| 52 | 2 - |M 2 M<?i/*D] 



r 2 [(r 2 ?r| 2 + Am 2 ] 131 " /ll|2(lff2 ' 2 " |/l2 ' 2) " r 2 [(r + rt™+ a™^ 1 " ^l 2 ^ 2 ) 

^r 2 + Z 2l )^:r^Am 2 ] ^ 2 + fe2 ' 2 ^' 2 -'^' 2 ) 
r(r2 + r) - W r |»-fc|'(M'-N') 



(r 2 + Am 2 )[(r 2 + r) 2 + Am 2 ]' 



2Am(r! + 2r) , ... 

+ (r 2 + Am 2 )[(r 1 + r) 2 + A TO 2 ] |g2 + " 2 " 9 <**i) 

2Am(r 2 + 2r) , 2 

" (r 2 + Am 2 )[(r 2 + r) 2 + Am 2 ] L92 - H A9lhl} 

4TAm 

^(gihig 2 h 2 ), 



(r 2 + Am 2 ) 2 

2 = 2fs W^l 2 + \ h i\ 2 M92K) + (M 2 + \h 2 \ 2 )Q( gi hl) + 2Q( gi h* l92 h* 2 )] 

+ ^1 [(l5i| 2 + Ih^Mg.K) + (\g 2 \ 2 + \h 2 \ 2 )-s( 9l hl) - 2% 9l h\g 2 h* 2 )\ 

r 2 - a™ 2 

+ (p + Am2)2 [(lgi| 2 + \hi\ 2 m 92 h* 2 ) + (\g 2 \ 2 + IhtftSbnht)] 
- ^P + Z^^^Am 2 ] ^ 2 ^ 2 ' 2 ^^) 

r(r 2 + r) - Am 2 2 
- 2 (r 2 + Am 2 )[(r 2 + r) 2 + Am 2 ] |g2 - /l21 
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Am 



T 1 [(T 1 +T) 2 + At 



r lsi + WM 2 -M 2 ) 



Af " rlffi-^i| 2 (l3 2 | 2 -|^| 2 ) 



r 2 [(r 2 + r) 2 + Am 2 ] 
r l[(ri 2( + r r)^ m ^ 1 + fel ' 25 ^^ 

^ r(k3 2 | 2 -|^i^| 2 ), 



(r 2 + Am 2 ) 2 

As = ^[l^l 2 - I3152I 2 - (M 2 - M 2 M 5 2^) - (I52I 2 - \h 2 \ 2 )3t( gi hl)} 
+ ^[l^l 2 - I51.92I 2 + (M 2 - l^il 2 )^^) + (M 2 - \h 2 \ 2 m gi h\)} 

or Am 

- "'g 1 \ 2 + \h 1 \ 2 )^{g 2 h* 2 ) + {\g 2 \ 2 + \h 2 \ 2 )^{g 1 h\)\ 



(r 2 + Am 2 ) 2 
T 2 - Am 2 



(r 2 + Am 2 ) 



(| ff ig 2 | 2 - \hM' 



Aa = ^r[(\9i\ 2 + N 2 )3(. 9 2^) + (M* + \h 2 \ 2 M 9l hl) + 2<Z( gi h* ig2 h* 2 )] 
Zl 1 

- ^[(l5i| 2 + N 2 )»(.92^) + (\g 2 \ 2 + |/i 2 | 2 )3(<7iM) - 2^{ gi h\g 2 h* 2 )\ 
2VAm 

' ~ r "g 1 \ 2 -\h 1 \ 2 Mg 2 h* 2 ) + (\g 2 \ 2 ^\h 2 \ 2 Mg 1 hl)} 



(r 2 + Am 2 ) 
T 2 - Am 2 

2 (rwF 3(9l ^ 2ft;) ' 



B " = sTTT 1 - 91 + ^ l|2 '- 92 " /l2 ' 2 + srk 1 - 91 " /l1 ' 2 '- 92 + ^ 2 ' 2 

- ^TW) [ (k|2 _ l^l 2 )^ 2 ! 2 - M*) + 43(51^)3(32^)] 

Bi = ^[l.9i^| 2 - |52^i| 2 + (kl 2 - M^G^S) + (l.92| 2 - \h 2 \ 2 n{ gi h\)} 

- ^[I.9i^ 2 l 2 - \92hi\ 2 (|.9i| 2 - N 2 )KGfc/iS) - (I.92I 2 - M 2 )»($i^)] 

- p + Am 2 [ (l - 91 ' 2 ~ l^l 2 )^ 2 ^) + d5 2 l 2 - M'Msi^)] - r^Am 2 ) ^ 1 ^ 92 ^' 
S2 = [290/1^5^2) + (Iffil 2 + M^M) + (I.92I 2 + \h 2 \ 2 m gi h\)\ 

~ 2lk [ " 2 ^iKg 2 h 2 ) + (1/nl 2 + |. 9l | 2 )3(.92^) + (\h 2 \ 2 + \g 2 \ 2 M gi hl)} 



T 2 + Am 2 
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[(| 5l | 2 + |/H| 2 )9(52fcS) + (|.g 2 | 2 + \h 2 \ 2 )-s( gi hl)} 



r(r 2 + Am 2 ) 
1 

2rr 2 
1 



fi 3 = ^[l5ig 2 | 2 - l^i^l 2 + (Iffil 2 - N 2 Mg 2 /* 2 ) - (M 2 - \h 2 \ 2 Mg 2 h* 2 )} 



nTi [l^ 2 | 2 - |.g l52 | 2 + (|. 9l | 2 - M 2 )H(«fe^) - (|.g 2 | 2 - N^M)] 

- r(r2 ^ m2) [(M 2 + \hi\ 2 Mg 2 h* 2 ) - (|.g 2 | 2 + \h 2 \ 2 M 9l hl)] 

+ r 2 + 1 Am2 (lgi-9 2 l 2 -l fe ^ 2 l 2 )^ 

B 4 = 2ffr [29( 5 i/ii02/»S) + C I 1 2 + N 2 )3(<7i^) - (k| 2 + N 2 )S*Gfc^)] 

- ^-[2^{ gi h\g 2 h* 2 ) + (|. 9l | 2 + |^| 2 )3(. 92 ^) - (| ff2 | 2 + \h 2 \ 2 m 9l h\)} 

- r(r2 ^ m2) [(kl 2 - M 2 Ms^ 2 ) - (l.g 2 | 2 - N 2 M.9i^)] 

Co — Ao,Cs — A3 , C4 — A4, 

c i = aff [M2I 2 - | 5 i5 2 | 2 - (I51I 2 - N 2 )R(<fefc$) - (|g 2 | 2 - \h 2 \ 2 m 9l hl)} 
2f| [M2I 2 - lsis 2 | 2 + (I51I 2 - \hi\ 2 W{g 2 h* 2 ) + (|g 2 | 2 - N 2 M<?i^)] 

T 2 — AT77 2 

■ [(M 2 - N 2 )»Gfcfc2) + (M 2 - N 2 )^/^)] 



2 



(r 2 + Am 2 ) 

+ ^[(r.+Vy+Am 2 ] 1 - 92 + /l2|2(l51 ' 2 - " l|2) + rfT^ te + " 2|2 ^ l/lI) 
T '-' " 1 | 52 - M 2 (M 2 - M 2 ) - F r ^ ,1 A T, a 2 J g 2 ~ h 2 \^( 9l hi) 



r 2 [(r 2 + r) 2 + Am 2 ] iaz Zl < 11 ' r 2 [(r + r 2 ) 2 + Am 2 ] 
2 

+ 

r(r 2 + r) - Am 2 



r(ri + r) - Am 2 , l2/l l2 ., l2 . 

+ ( r 2 + Am 2 )[ ( r 1 + r) 2 + Am 2 ] |gl + - 1/121 } 



(r 2 + Am 2 )[(r 2 + r) 2 + Ar 



-Isi-MW-M' 



2Am(r 2 + 2r) , 2 

- (r 2 + Am 2 )[ ( r 2 + r) 2 + Am 2 ] |gl - ^ 5(52 " 2) 

2Am(r! + 2r) ., , 

+ (r 2 + Am 2 )[(r 1 + r) 2 + Am 2 ] |gl + h ^ ^ 92h2) 

- (r2 + Am2)2 ^(giM 2 fe 2 ) ; 
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C 2 = ^2 [d.9i| 2 + \hi\ 2 )%g 2 h* 2 ) + (|.g 2 | 2 + \h 2 \ 2 )Z(g 1 h* 1 ) + 2Z(g 1 h* 1 g 2 h* 2 )] 
+ tL [d.9i| 2 + \hi\ 2 m92h* 2 ) + (\g 2 \ 2 + \h 2 \ 2 m gi hl) 2Z( 9l h* l92 h* 2 )] 



2Ti. 



(r r 2 2 + ^ 2 2 )2 [(l.9i| 2 + \hA 2 )%g 2 K) + (\g 2 \ 2 + \h 2 \ 2 M 9l hl) 



1.92 + h 2 \ 2 Q( 9l ht) - - =- \ i' ~ 21 |g 2 - /i 2 | 2 3(MD 



ri[(r! + r) 2 + Am 2 ] |yz 1 Z| vyi i; r 2 [(r 2 + r) 2 + Am 2 ] 

Am(Fl + ^ r|* + M 2 (| ff2 | 2 - N 2 ) - 7T^W(| 5152 | 2 - I/mM 2 ) 



(r 2 + Am 2 )[(r! + r) 2 + Am 2 ] 131 11 1 Z| ; (r 2 + Am 2 ) 

Am(r 2 + 2r) , 2 l2 l2 , 

+ (r 2 + Am 2 )[(r 2 + r) 2 + Am 2 ] |gl - /l11 (l52 ' - 1/121 } 

r(r 2 + r) - Am 2 2 
- 2 (r 2 + Am 2 )[(r 2 + r) 2 + Am 2 ] |gl - h ^ ^ 92h2) 

rb-^ 2 | 2 (| 5 i| 2 -|^| 2 ) 



r 2 [(r 2 + r) 2 + Am 2 ]' 
+ r 1 [(r 1 + rF + A^] l92 + fe2|2(lfll|2 - |fel|2) ' 

Do = ^r\9i - ^i| 2 | 52 + h 2 \ 2 + g^r\9i + M 2 |S2 - h 2 \ 2 

[(I51I 2 - I^il 2 )(l5 2 | 2 - \h 2 \ 2 ) + mgihlMg2hZ)] 



4(r 2 + Am 2 ) 



Am ■ [(| ffl | 2 - N^M) - (\g 2 \ 2 - \h 2 \ 2 M 9l hl)] . 



2r(r 2 + Am 2 ) 



Di = ^[l3i^2| 2 - Ig^l 2 + (| 5l | 2 - Ih^Mg.h*) + (\g 2 \ 2 - \h 2 \ 2 M 9l hl)} 
1 [\gih 2 \ 2 - Ig.h,] 2 - (| 5l | 2 - Ih^Mg,^) - (\g 2 \ 2 - l^ 2 )^^)] 



2rr 2 
1 



[(| ffl | 2 - |/n| 2 M 52 /^) + (|.g 2 | 2 - \h 2 \ 2 m 9l hX)} 



T 2 + Am 2 

- r ( rw) 3(gift ' g2 ^- 

D 2 = ^[2%{ 9l h\ 9 * 2 h 2 ) + (|^| 2 + \ 9l \ 2 mg 2 K) + {\h 2 \ 2 + |g 2 | 2 )3( 5l ^)] 
1 [2%{ gi hlg* 2 h 2 ) {\h x \ 2 + \ 9l \ 2 Mg 2 h;) (\h 2 \ 2 + \g 2 \ 2 )-s( 9l hl)] 



2TT 



r 2 + Am 2 [(|gl12 + \ h ^92h* 2 ) + (\g 2 \ 2 + M 2 )^!)] 



Am -(k^2| 2 -| 5 2^i| 2 ), 



r(r 2 + Am 2 ) 
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D* 



1 



151521 



l/MM 2 + (M 2 - \hx\ 2 Mg 2 h* 2 ) - (\g 2 \ 2 - \h 2 \ 2 M 9l hl)] 



2TT 1 

' Wl 2 - \h x h 2 \ 2 - (|. 9l | 2 - M 2 )H(ff2ft$) + (| 52 | 2 - \h 2 \ 2 M 9l ht)] 



2rr 2 



A? 



[(| 5l | 2 + l^l 2 )^^) - (| 52 | 2 + l^l 2 )^!^)] 



D 4 = 



r(r 2 + Am 2 ) 

r2 A 2 {\9i92?-\hik 2 \% 

2 ^r[^(9ih*i9* 2 h 2 ) - + \gi\ 2 Mg 2 h* 2 ) + (\h 2 \ 2 + \g 2 \ 2 M gi hl)] 
1 



2rr 2 



[2%{ gi h\g* 2 h 2 ) + (l/jil 2 + \gi\ 2 Mg 2 h* 2 ) ~ {\h 2 f + | 52 | 2 )9( ffl ^)] 



Am 



r(r 2 + Am 2 ) 



[(| 5l | 2 - IfcHRG/afcS) - (| 52 | 2 - l^l 2 )^!^)] 
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